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Abstract 

The matrix model of Kapustin, Willett, and Yaakov is a powerful tool for exploring the 
properties of strongly interacting superconformal Chern-Simons theories in 2+1 dimensions. 
In this paper, we use this matrix model to study necklace quiver gauge theories with A/" = 3 
supersymmetry and U{Ny gauge groups in the limit of large A^. In its simplest application, 
the matrix model computes the free energy of the gauge theory on S^. The conjectured F- 
theorem states that this quantity should decrease under renormalization group flow. We show 
that for a simple class of such flows, the F-theorem holds for our necklace theories. We also 
provide a relationship between matrix model eigenvalue distributions and numbers of chiral 
operators that we conjecture holds more generally. Through the AdS/CFT correspondence, 
there is therefore a natural dual geometric interpretation of the matrix model saddle point 
in terms of volumes of 7-d tri-Sasaki Einstein spaces and some of their 5-d submanifolds. As 
a flnal bonus, our analysis gives us the partition function of the T{U{N)) theory on S^. 
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1 Introduction 



Exact results in strongly-interacting field theories are generally rare. In supersymmetric field 
theories, supersymmetry places strong constraints on various properties of chiral operators, 
and exact results pertaining to these operators might be possible even at strong coupling. 
For three-dimensional superconformal theories, recent progress in finding such exact results 
that hold at any coupling was made in (l}]3] , where the partition function of superconformal 
theories on with M >2 supersymmetry, as well as the expectation values of certain BPS 
Wilson loops, were reduced from path integrals to finite-dimensional multi-matrix integrals. 
This major simplification was achieved through the localization technique developed in [2] 
for four- dimensional theories. 

A consequence of this work is the realization that the "free energy" F defined as minus 
the logarithm of the path integral on S"^, 

F = -log|Zs3|, (1.1) 

with appropriate subtractions of power law divergences, might represent a good measure of 
the number of degrees of freedom in any field theory, supersymmetric or non-supersymmetric. 
One way in which F can be thought of as a measure of the effective number of degrees of 
freedom is the conjecture made in [5] that F decreases along renormalization group (RG) 
flows and is stationary at RG fixed points. This conjecture was called the "F-theorem" in |5] 
and was tested in a few RG flows in large supersymmetric U{N) gauge theories (5|[6]. If 
the F-theorem is true, then F would be a 3-d analog of the central charge c from two- 
dimensional field theory, which is known to have the same monotonicity property along RG 
flows |7], and it would resemble the Weyl anomaly coefficient a from 4-d theories, which is 
also believed to decrease along RG trajectories [8|. Actually, the free energy F also resembles 
a in another way: just like a, F can be used to find the exact R-symmetry in the infrared 
(IR) by computing F as a function of a set of trial i?-charges and then maximizing it [2|. 
The analogous procedure in 4-d theories is called "a-maximization" fo). 



In the context of the AdS/CFT correspondence 10 12 , the free energy F can also be 
computed holographically from the gravity side of the correspondence. In particular, for a 
CFT dual to AdS^ of radius L and effective four-dimensional Newton constant Gn, F is 
given by |13j 

F = — . (1.2) 
2Gn 
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It was shown in fH] that in any CFT F can also be interpreted as an entanglement entropy 
between a disk and its complement in the M^'^ theory. In turn, this entanglement entropy 



equals the holographic a* function defined in 15 , 16 that was shown to always decrease along 
holographic RG fiows. 

If the 74^5*4 background mentioned above arises as a Freund-Rubin compactification 
AdS^ X y of M-theory, where F is a seven- dimensional Sasaki-Einstein space threaded by 
units of four-form fiux, then the quantization of the AdS4 radius in Planck units implies 



that at large N eq. (1.2) becomes 17 



Here, the volume of Y is computed with an Einstein metric that satisfies the normalization 
condition Rmn = Qgmn- The Freund-Rubin solution AdS^ x Y arises as the near- horizon limit 
of a stack of N M2-branes placed at the tip of the Calabi-Yau cone X over Y. The Ar3/2 
behavior of the number of degrees of freedom had been known for quite some time, as the 
same large N dependence appears in other quantities such as the thermal free energy that was 



computed in 18 more than ten years ago. A field theory explanation of this peculiar large N 
dependence had been lacking until recently, mostly because explicit Lagrangian descriptions 
of the field theories living on coincident M2-branes have been discovered only in the past few 



years. Starting with ABJM theory 19 that describes A^ M2-branes sitting at an orbifold 
singularity of C^, there are now many Chern-Simons matter U{N) gauge theories that are 
proposed to describe the effective dynamics on A^ M2-branes placed at the tip of various 



Calabi-Yau cones (see for example 20 24]). However, only few of these dualities have been 
extensively tested. Extensive tests are difficult to perform because supergravity on 74^5*4 x Y 
is supposed to be a good approximation to the dynamics of the CS-matter gauge theories 
only as one takes the gauge group ranks to infinity while keeping the CS levels fixed. In this 
limit, the 't Hooft coupling N/k becomes large, and there are no perturbative computations 
that one can perform. 

That in A/" > 2 theories one can write F exactly in terms of a matrix integral means 
that by evaluating this integral one can test some of these AdS^/CFT^ dualities that have 
been put forth in recent years. In particular, one can provide a field theory derivation of the 
F oc N^/"^ large A^ dependence of the number of degrees of freedom (1.3) on A^ coincident 



M2-branes. Moreover, one can compare the coefficient of N^^'^ in (1.3) to the volume of the 



internal space Y that one can compute independently by integrating the square root of the 
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determinant of the Sasaki-Einstien metric on Y. Such comparisons were made in 25 in the 



case of ABJM theory, in 17 ,26 for a large class of theories with A/" = 3 supersymmetry, and 
in [5,27,28 for many theories with J\f = 2 supersymmetry. In 25,26 F was computed as a 
function of the 't Hooft coupling, while the other works focused only on the strong coupling 
regime and used the method developed in [17j to evaluate the multi-matrix integrals at large 
N and fixed CS levels in a saddle point approximation. 

In this paper we build upon the work in [5 , 17 in several ways. Our first main result 
consists of an infinite class of RG flows whose IR and ultraviolet (UV) fixed points preserve 
Af = 3 SUSY. Via the AdS/CFT correspondence, at large such a flow is dual to a 
holographic RG flow between two AdS^ x Y extrema of 11-d supergravity. That the extrema 
preserve A/" = 3 supersymmetry means that the cones over the spaces Y are hyperkahler 
and, because of that, the spaces Y are called tri-Sasakian. (We will provide a more detailed 
exposition of hyperkahler spaces at the beginning of section 2.1[ ) The 11-d SUGRA solutions 
that represent M2-branes sitting at the tip of a hyperkahler cone were constructed in [29j, 
where they were also related through string duality to a brane construction in type IIB string 
theory. The type IIB brane construction consists of A^ D3-branes filling the 0126 directions, 
and a sequence of {pa, Qa) five-branesj^ 1 < a < d, filling the 012 directions and sitting at 
fixed angles in the 37-, 48-, and 59-planes. The {p, q) branes break the D3-brane stack into 
d segments, and the three-dimensional U{NY CS-matter gauge theories then live on the 
segments. The A/" = 3 RG flows we want to study correspond to removing one of the (p, q) 
branes or to a brane and a {p2,q2) brane combining into a {pi + P2,qi + ^2) brane. 

We find that the F-theorem is satisfied in these examples. 

Our results for the F-theorem follow most naturally from the gravity side of the AdS/CFT 
correspondence which provides an SL{2,'R) invariant result for F as a function of arbitrary 
{pa, Qa)- We provide field theory confirmation of the gravity result in the cases where Pa = 
or Pa = 1 for each a. While the 3-d field theory is simplest when = or 1, for pa > 1, 
ref. |30] provides a more complicated field theory description involving an interpolating 
T{U{N)) theory that implements a sort of local 5* duality. Instead of using the T{U{N)) 
theory, we derive the matrix model for > 1 by bootstrapping from our large A^ results. 
We check that the matrix model yields the correct answer in the large A^ limit, and also that 
it is invariant under S'L(2,Z). As a bonus, we discover the matrix model of the T{U{N)) 
theory]^ The matrix model for arbitrary pa is our second main result. 



^We adopt the convention where by a (p, q) five-brane we mean a brane with p units of NS5 charge and 

q units of D5 charge. 

^As we were completing this work, 31 appeared which contains this same result. Also, M. Yamazaki, 
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Along the way we tie several loose ends left off from 17 . In 17 the matrix model was 
solved explicitly only for theories on (1, g^) branes with a < 4. We provide a general solution 
that holds for any number of {paiQa) branes. In |17| it was checked numerically in a few 



examples that the M-theory prediction for the volume of Y that can be extracted from (1.3) 



agrees with the geometric computation performed by Yee 33 . We build on Yee's work and 
prove that in general the two computations of Vol(y) agree. In 17 it was conjectured that 
the volume Vol(F) can be expressed in terms of a certain sum over trees. We provide a proof 
of this tree formula. 

Our third main result, contained in Section [3} is a relationship between the eigenvalues 
in the matrix model and the number of chiral operators in the field theory. Define il){r) to be 
the number of chiral operators with i?-charge smaller than r for the = 1 gauge theory. The 



authors of 34,35 demonstrated that there is a relationship between Tp{r) and the volume of 



the Sasaki-Einstein space in the large r limit. Given (1.3), there must also be a relationship 
between il){r) and F. In fact, as we show in this paper for the necklace theories, more precise 
relationships can be established between the matrix model and operator counting problems. 
The operators in the necklace theories also have a monopole charge m corresponding to m 
flux units in a diagonal subgroup of the gauge group. Thus, we may consider ijjir^m) to be 
the number of operators with i?-charge less than r and monopole charge less than m. Let 
Xab be a hypermultiplet transforming under the fundamental of the 6th gauge group and the 
antifundamental of the ath gauge group. We can consider ipXabi''^^ ^) ^'^ t)e defined as above 
but now with the operator Xab = 0. In the large N limit, the matrix model is solved by a 
saddle point approximation for which the eigenvalues are complex numbers Aq = N^^'^x + ii/a- 
The eigenvalues can be parametrized by an eigenvalue density p{x), which turns out to be 
the same for each gauge group, and an imaginary part ya{x). Our two results are 



dr'^dm 



drdm 



= -Pix), (1.4) 

m=rx/ii P 

r 

= -p{x)[yb{x) - yaix) + R{Xab)] , (1.5) 

m=rx/ii P 



where p = SF/AtcN^/"^. (H ere we take the liberty of replacing the operator counts, which 
are discrete functions, with continuous approximations.) We believe these relations will hold 
more generally. 



T. Nishioka and Y. Tachikawa have informed us they have independently derived the T{U{N)) partition 



function 32 
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2 Volumes of Toric tri-Sasaki Einstein Spaces 



On the gravity side of the AdS/CFT duahty, we have M-theory backgrounds generated by 
placing a stack of M2-branes at the tip of a hyperkahler cone, with large. In this section 
our aim is to compute the free energy of the M2-brane theory purely from the supergravity 



side of the correspondence. We start by introducing in section |2.1| the main ingredients 
in constructing the 11-d supergravity solution, namely the toric hyperkahler spaces. In 
section 2.2 we build upon the results of Yee 33 and express the volume of these spaces in 



terms of the volume of a certain polygon for which we will provide a field theory interpretation 



later on. In section 2^ we comment on some field theory implications of this formula, and 
show explicitly that the free energy decreases along certain RG flows, in agreement with the 
F-theorem proposed in [s]. 

2.1 Toric Hyperkahler Cones from a Quotient Construction 

We start by introducing the toric hyperkahler cones. The following discussion draws heavily 



from 36 and 33 




A hyperkahler manifold possesses 4n real dimensions and has a Riemannian metric g 
which is kahler with respect to 3 anti-commuting complex structures Ji, J2, and J3. These 
Ji furthermore satisfy the quaternionic relations = J| = J| = J1J2J3 = —1- The 
simplest example of a hyperkahler manifold is the four-dimensional space of quaternions 
M = M^, endowed with the standard line element. A single quaternion g G EI = can be 
represented two-by-two complex matrix 



(2.1) 



parametrized in terms of two complex variables u and v. In terms of u and v the metric is 
ds'^ = \duf + \dvf and the three kahler forms are 

= --(duAdu + dvAdv), (2.2) 
{u^ - iu^) = i{du A dv) . (2.3) 

The u'^ transform as a triplet under the SU{2) symmetry that acts as left multiplication on 
The quaternions EI can also be written as a U{1) bundle over where the U{1) fiber 
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shrinks to zero size at the origin of M^. This description comes from the uphft of the Hopf 
fibration from S*^ to and makes exphcit an 5*0(3) x U{1) subgroup of the 0(4) rotational 
symmetry of W^. More exphcitly, if one writes 

M = cos (6/2) e-^-^/^+i^ , = sin (6/2) e'^'^/^-iV' ^ (2.4) 
then the standard hne element on EI becomes 



ds =dudu + dvdv = - \- r [ dip cos(^)c?0 , (2.5) 

4 r \ 2 J 

where r = (r sin 6 cos 0, r sin 6 sin 0, r cos 6) is the usual parameterization of in terms of 
spherical coordinates. The coordinate ip G [0, 2tt) parameterizes the U{1) Hopf fiber. From 



eq. (2.1) we see that rotations in the U{1) fiber correspond to phase rotations of q. 

Another example of a hyperkahler manifold is the Cartesian product of ci-copies of the 
quaternions M.'^ = M^'^, also considered with the flat metric 



a=l 



Idr^ ( 1 ^ ^ 

1^ + ( #a - 2 COs{9a)d(f)a 



(2.6) 



Starting with M.'^, one can construct construct a large number of further examples of toric hy- 
perkahler spaces using the hyperkahler quotient procedure of Hitchin, Karlhede, Lindstrom, 



and Rocek 37 . A toric hyperkahler manifold can be defined to be a hyperkahler quotient 



of EI for some integer dhy a {d — ?T,)-dimensional subtorus C T . The underlying reason 



that the quotient continues to be hyperkahler is that the triplet of Kahler forms (2.2) is 
invariant under the U{1) action that sends u — ?■ we*" and v — )■ fe~*". 

The hyperkahler quotient procedure requires the data of how A^ = T*^-" sits inside T'^. 
The inclusion N G T'^ can be described by the short exact sequence of tori 

— y T'^-" ^ yd ^ — ^ g ^ (2.7) 

where we also introduced the quotient = jrpd-n^ ^^jj ^-^^ inclusion map i and the 
projection map vr. Each torus in this sequence can be identified with IR'^/(27rZ'^), so the 
data of how A^ sits inside T'^ can be encoded in how the standard lattice Z*^"" C W^~^ sits 
inside 7/ C M'^, or equivalently how 7/ C projects down to Z" C M". Such a construction 
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can be described by a short exact sequence of vector spaces 



nd—n ^ , Tad ^ 



which restricts to a short exact sequence of free-Z modules (lattices): 



— y Z^"" ^ ^ Z" — ^ . 



{21 



(2.9) 



Since Q and /3 are linear maps, they can be represented by matrices: Q by a c? x [d — n) 



matrix and /3 by an n x c/ matrix. That the sequence (2.8) is exact means that precisely 
n columns of (3 are linearly independent (i.e. (3 is surjective), that (3Q = 0, and that the 



d — n columns of Q are linearly independent (i.e. Q is injective). That (2.8) restricts to ( |2.9 ) 
further implies that Q and /3 have integer entries. The torus N = M'^^"/ (27rZ'^^") as we 
defined it has volume {2ttY^"'. 

The hyperkahler quotient of M.'^ by the torus is defined to be the zero locus of a set of 
moment maps in addition to a quotient by the torus action. The 3{d — n) moment maps are 
compactly expressed using the matrix equation 



a=l 
d 

= E« 



a=l 



ra3 - ira2 



+ 



A, 



a'i 



\al - iX 



a2 



Tal + ^ra2 



Aal + «A, 



a2 



-A, 



a3 



(2.10) 
(2.11) 



where cxi, (T2, and are the Pauli spin matrices. The hyperkahler quotient X is then 

X = W^///N = fi-\0)/N . 



(2.12) 



We are particularly interested in the case where Aa = and the corresponding hyperkahler 
manifold is a cone. The base of this cone is a 4r;, — 1 real dimensional Riemannian manifold 
with positive scalar curvature and a locally free action of SU{2) that descends from the 
SU{2) rotating the three complex structures. The base of such a cone is called a tri-Sasaki 
manifold. Indeed, a Riemannian manifold {Y, g) is tri-Sasaki if and only if the Riemannian 
cone X = (]R+ x Y, dr"^ + gijdy^dy^) is hyperkahler. The induced metric on X from is 
Ricci flat, which in turn implies that the induced metric on Y satisfies Einstein's equations 
with a positive cosmological constant. In other words, Y is also an Einstein manifold. 
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2.2 The Volume of the tri-Sasaki Einstein Base 



Given a toric hyperkiihler cone X, we would like to compute the volume of the base Y with 
respect to the induced metric from M.'^. A brute force approach would be to introduce the 



variables tj such that 



5^/3a/i. (2.13) 



ra 

These constrained automatically satisfy the moment map conditions. Plugging this ex- 



pression into the line element (2.6) on yields a metric on X x T*^ Integrating the 
square root of the determinant of the metric up to a finite radial coordinate 1 = and 
dividing by VoKT"^"") would yield Vol(F)/4n. 

The approach outlined above appears to be difficult. Instead, we will use a result due to 
Yee [33], which was proven using an elegant localization argument: 

Yee's Formula. Consider a tri-Sasaki Einstein manifold Y defined via the short exact 



sequence (2.9) and hyperkdhler quotient construction described above. If the metric is nor- 



malized such that Rij = 2{2n — l)gij, then the volume is 

od— n+l _2n 

Vol(F) = ..,..„w.n I \\\d<Pn\ I — r^. (2.14) 




(2n - 1)! Vol (N) J ^ J-1 / 



This volume formula has the forgiving property of being invariant under rescaling the 
matrix Q — )■ XQ. Under this rescaling the torus volume changes, Vol(A^) — )■ A"~'^Vol(A^). 
The factor A"^*^ is then canceled by the Jacobian introduced upon rescaling — 0/A. We 
discussed above that the columns of Q can be chosen such they form a Z-basis of the kernel 
of P and Vol(A^) = (27r)'^~". Note that if we were not so clever in our choice of Q, the volume 
formula would still give us the correct answer because of this scaling invariance. 

So far we have assumed that N is isomorphic to T'^, but we can also choose N to be 
isomorphic to Nt x T*^ for some finite abelian group Nt- In that case, the image of /3 is no 
longer all of Z". The cokernel of /3 is dual to Nt- The volume of N for our choice of Q is 
therefore (27r)'^-" \Nt\ = (271)^^"" |Coker/3|. 

We now rewrite this volume integral in a more convenient fashion. Let us assume that we 
chose the columns Q to form a Z basis of Ker(/3). In this case, the volume formula simplifies: 



.3n—d 




Vo.(V)^^/MM^|M^^^^. (2.15) 
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A well-known result from Fourier analysis is 



1 + 



dye 



-\y\+ixy 



(2.16) 



Then Vol(y) can be rewritten as 



Vol(r) 



TT' 



3n—d 



{2n - 1)! ■ 2^-1 




He-^y-^dy, exp J^^Q^^V • (2.17) 



a=l 



k,a 



We are allowed to switch the order of integration, integrating over the (jP first We obtain 



Vol(r) 



TT 



2n 



d~n 



(2n- 1)! -2 



rt— 1 



(2.18) 



Ka=l 



k=l 



To integrate over the delta functions, note that we can get a basis for 7/ by taking the 
basis for Ker(/3) and pullbacks of the basis for Im(/3). The Jacobian for transforming from 
the standard basis of to this basis must be one, since both bases generate the same 
lattice. In our new coordinates, we have two kinds of variables: Si corresponding to the 
columns of Q and ti corresponding to the rows of (5. The product of delta functions is just 
5{si)5{s2) ■ ■ ■ 5{S(i-n) and can be performed straightforwardly. If the rows of (5 span Z", the 



tj can be written 



The integral reduces to the following useful form: 



(2.19) 



Voi(r) 



TT 



2n 



(2n - 1)! ■ 2"-i 




exp 



E 

a=l 



(2.20) 



Note that if the rows of /3 do not span Z", then the Jacobian will have an extra factor of 



I Coker(/3)|, which cancels the | Coker(/3)| in Vol(A^). So (2.20) holds even if the rows of (3 
do not span Z". This integral form provides us with a corollary to Yee's Formula: 



The integral (2.17) is not absolutely convergent. 



However, if we multiply the inte- 
grand by exp [— efc(0''")^] , for some small ej; > 0, then the integral does converge abso- 
lutely. We can then take the limit 0+. If we integrate out the Ua then we get 



lim,^^o+ ji 



2? 

(2' 



n-l)\ J Ulj = l 



1+(E..Q^0'=) 

integrate out the (j)^ then we get lim,^^o+ (2«-T)!^2"-i / (na=i 11^=1 
J2^,^S^ I (nti e-\y^\dya) nti S (Ea QlVa). 



C^P-{l2a QkVci) /4gfc _ 
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Corollary 1. i) If a column Pa is removed from /3, Vol(y) increases, ii) If two columns 
(3a and f3b of (5 are combined to form the new column f3a + Pb of a new f3' with one fewer 
columns, Vol{Y) either increases or stays the same. The volume remains the same if and 
only if the two columns are proportional, Pa = cPb for some c G . 



Proof, (i) is true because the integrand of (2.20) (which is positive definite) increases when 
a column of P is removed, (ii) is true because the absolute value of a sum is less than or 
equal to the sum of the absolute values. Equality occurs only when the expressions inside 
the absolute values are proportional. □ 



The volume formula (2.20) leads to another useful corollary of Yee's Formula: 



Corollary 2. Let P G be the polytope 



V 



t G 



a=l 



< 1 



(2.21) 



Let S"*" ^ be a {An — 1)- dimensional sphere with unit radius. The volume of the tri-Sasaki 
Einstein manifold satisfies the relation 



Vol(y) 
Vol(54"-i) 



Vol(P) . 



(2.22) 



Proof. Introduce the notation a = Ylt=i X]j=i f^ajt 



with a > if at least one of the tu is 



non-zero, and write ti, = ahu for some new variables hh. The hk are constrained to live on 



the boundary of V. Eq. (2.20) can be rewritten as 



Vol(F) 



vr 



2n 



{2n - 1)! ■ 2^ 



\k=l 



va=l 



1 I j daa'^-^e-''. 



The volume of V can be also written in terms of these variables: 

d 



Vol(P) 



,a=l 



5Z ^-J'^J' 



1 I j daa""-^ 



Performing the integrals in a and comparing the last two formulas, we obtain 



Vol(y) 



{2n - 1)! ■ 2'- 



-Vol(P) 



(2.23) 



(2.24) 



(2.25) 
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or, using the fact that the volume of the (4n — l)-sphere is Vo\{S^^ ^) = 27r^"/ {2n — 1)! the 
desired result follows. □ 

A final corollary is an explicit result for the volume of the tri-Sasaki Einstein spaces 
relevant for the gauge theories we discuss below: 

Corollary 3. In the case n = 2, choose /3 such that the two-vectors Pa He in the upper half 
plane and order them such that /3a A /3a+i > Oj^ The volume of the tri-Sasaki Einstein space 
Y IS 



d 



Vol(F) = -y ^^^^ (2.26) 

where we have defined the quantities ■jab = IPa /\ Pb\, o-a = Ylb=i '^a.b, and f3d+i = —Pi- 

Note that the volume of Y is independent of the sign of the Pa and their order inside 
p. The order of the Pa in the corollary is their order around MP^. (Flipping the sign of 
any given column of P does not change the hyperkahler quotient M./ / /N because the torus 



defined by (2.9) is invariant under such a sign flip. Therefore it is natural to identify Pa 
with —Pa and consider the ordering of the Pa = {Pa, Qa)^ to be the ordering of \pa, Qa] G MP^ 
around MP^) 

Proof of Corollary We split the integral into various regions according to whether Pajtj 
is positive or negative. There are 2d such regions. Since the integral does not change if 
we replace each tk with —tk, we only need to consider d regions. The ordering of the Pa 
guarantees that any two consecutive columns determine a region boundary. 

Now we choose a and compute the integral in a region bounded by the lines u = 
— J2jt^ajtj = and V = P(a+i)jtj = 0. The prefactor in the integral is 7r^/12, but 
we need to multiply by two since there are two such regions. We get 



l-OO l-OO / P 

/ du dv exp — y 

6 1a(a+l) Jo Jo \ b^ 7a{a+l) 



TT^ 7a(a+l) 



(2.27) 



6 (EbTafe) (Efe7(a+l)fe) 

Summing the regions yields the volume stated in the corollary. □ 
Given Corollary |2| there exists an equivalent proof that involves computing the area of 



V from the definition (2.21) 



^We define (?) A r, 1 = ad - 5c. 
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{P2, q2 

(pi, gi) 5-brane 



{Pd, qd) 




N D3's 



Figure 1: A schematic picture of the brane construction. The N D3-branes span the 0126 
direction, and the {paAa) 5-branes span the 012 directions as well as the lines in the 37, 48, 
and 59 planes that make angles 6a = arg(pa + iqa) with the 3, 4, and 5 axes, respectively. 
The three-dimensional A/" = 3 theories considered in this paper live on the 012 intersection 
of these branes. 



The authors of 17 conjectured a formula for the volume of these n = 2 tri-Sasaki Einstein 



spaces. Their formula is interesting because it does not rely on an ordering of the Pa and 
makes the permutation symmetry of the columns of /3 manifest. We have been able to 
promote this conjecture to a theorem. As the techniques for the proof are not typical of the 
main arguments in the paper, we include the proof as Appendix [Aj 



Tree Formula. The area of the polytope V described in (2.21) in the case n = 2 can be 
written 

X^(V",E)eT n(a,f))e-E 7a6 



Area(P) = 2- 



(2.28) 



Ua=l 

where T is the set of all trees (acyclic connected graphs) with nodes V = {1,2, .. . ,p} and 
edges E = {(ai, hi), (02, 62), • • • , (op-i, VOI- 



2.3 Brane Constructions and an F-theorem 

Consider the following brane construction in type IIB string theory. A stack of coincident 
D3-branes spans the 01236 directions with the 6 direction periodically identified. Let there 
be bound states of NS5- and D5-branes. Denote the number of NS5- and D5-branes in 
the bound state by pa and qa respectively. These {pa, ga)-branes intersect the D3-branes at 
intervals around the circle and span the 012 directions. Each {pa, ga)-brane lies at an angle 
9a in the 37, 48, and 59 planes where 9a = arg(pa + iqa). These brane constructions are 
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known to preserve 6 of tlie 32 supersymmetries of type IIB string theory 38 ,39 ^ 

This brane construction cannot be described rehably within type IIB supergravity be- 
cause the dilaton becomes large. A better description can be obtained after a T-duahty along 
the 6 direction and a lift to M-theory, where the resulting configuration can be described 
within 11-d supergravity. The geometry depends on A^. When is small, the geometry is 



d2,1 



X X where X is the hyperkahler cone discussed above for n = 2 29 . In the large A^ 
limit, the D3-branes produce a significant back- reaction, and close to them the geometry is 
AdSi X Y where y is a tri-Sasaki Einstein space 20 . In both cases, the charges {pa, qa) = 0^ 
are the columns of /3. 

In the case where pa = or 1, this brane construction admits a simple A/" = 3 super- 
symmetric 2+1 dimensional field theory interpretation. The (1, ga)-branes, a = 1, . . . ,d, 
break the D3-branes up into d segments along the circle. For each segment, we have a U{N) 
Chern-Simons theory at level ka = qa+i — Qa- The (0, gj)-branes, i = 1, . . . , np, also intersect 
the D3-branes. At each intersection, we have massless strings that join the D3-branes to 
the (0, gi)-brane that correspond to flavor fields in the fundamental representation. (See 



figure 2 ) These gauge theories are described in greater detail in 20,40 



As described in the introduction, it has been conjectured that the logarithm F of the 
partition function of the Euclidean field theory on S*^ serves as a measure of the number 
of degrees of freedom in the theory, being an analog of the conformal anomaly coefficient a 
of a 3+1 dimensional field theory. While the a-theorem is a conjecture that the conformal 
anomaly a decreases along RG flows, the conjectured F-theorem states that F should de- 



crease along RG flows [5]. As given in (1.3), the AdS/CFT correspondence predicts that in 
the large A^ limit, F ~ l/^Vo\{Y). Thus, along RG flows, Vol(y) should mcrease. 

It is remarkable that for the class of theories we consider here. Corollary [T] confirms this 
expectation that Vol(y) should increase along RG flows. The UV of the field theory should 
correspond to looking at our D3-brane construction from a great distance as only high energy 
excitations will be able to get far from the branes, while the IR should correspond to getting 
very close to the D3-branes. The simplest realization of such an RG flow is to add a mass 
term for flavor fields corresponding to a (0, gj)-brane. Such a mass term corresponds to 
introducing a small distance between the (0, gj)-brane and the stack of D3-branes. In the 
UV, these flavor fields will contribute to F, while in the IR, the flavors and corresponding 
(0, gj)-brane should be absent. 

A slightly more complicated realization involves an intersecting (0, g)-brane and a (1, 0)- 

^Note that a brane with charges {—Pa, —<la) is an anti-(pg, (7a)-brane rotated by 180 degrees, which is the 
same as a {pa, ga)-brane. The [pa, qa) charge is most naturaUy defined up to this overaU sign. 
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brane. We have a U{N) vector multiplet on each side of the (1, 0)-brane. The (0, g)-brane 
produces q fundamental flavor fields charged under one U{N) and q anti-fundamentals 
charged under the other. Adding real masses of the same sign for each flavor field corre- 



sponds to a "web deformation" 39 where close to the D3-branes the (0, q)- and (1, 0)-brane 
form a (1, g)-brane bound state. Far from the D3-branes, because we have fixed boundary 
conditions for the (0, g)- and (l,0)-brane by specifying their angles, the two branes will re- 
main separate. From a large distance, we will not see that the two 5-branes have combined. 
We will only see their asymptotic regions where they remain separated. Closer up, we will 
see that the branes have made a bound state. From the point of view of our volume formula, 
combining these two branes increases the volume and thus decreases F. 

In fact. Corollary [l] and our volume formula put no restriction on the type of {p, g)-brane 
we remove or the type of {p, q)- and {p', g')-brane we combine to form a {p + p',q + q') bound 
state. In all cases, the volume will increase, corresponding to a decrease in F. Unfortunately, 
from the field theory perspective it is not clear in general to what these more general types 
of RG flow correspond. 



3 Field Theory Computation of the Free Energy 
3.1 A/" = 3 Matrix Model 

Consider the Af = 3 supersymmetric 2+1 dimensional Chern-Simons theories corresponding 
to the {p, g)-brane constructions described above where p = or 1. Let there be d U{N) 
gauge groups at level ka = qa+i—qa, matter fields Aa and Ba in conjugate bifundamental rep- 
resentations of the (a — l)st and ath gauge group and Ua pairs of flavor fields transforming in 
fundamental and anti-fundamental representations of the ath gauge group, with a = 1, . . . ,d 
(see figure |2|. As explained in [l], the partition function for these necklace quivers localizes 
on configurations where the scalars cTq in the vector multiplets are constant Hermitian ma- 
trices. Denoting the eigenvalues of cTq by A^^j, I < i < N, the partition function takes the 
form of the eigenvalue integral 

Z = j (^dXa,}j L,{{XaA)Lmi{Ki}) (3.1) 
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Figure 2: A necklace quiver gauge theory where the gauge sector consists of d U (N) gauge 
groups with Chern-Simons coefficients ka- The matter content consists of the bifundamental 
fields Aa and Ba, as well as pairs of fundamentals and anti-fundamentals transforming 
under the ath gauge group. 



where the vector multiplets contribute 



= M n ( n 2 sinh[7r(A„,i - Xaj)] j exp | in ^ kaXli j 



(3.2) 



and the bifundamental and fundamental matter fields contribute 

2cosh7rAa, 



n n 

a=l \ i,j 



2cosh[7r(Aa-i,i - Xa,j)] 



(3.3) 



We follow the recipe suggested in |17j for analyzing this matrix model in the large N 
limit. We write A = N^^'^x + iy and assume that the density of eigenvalues p{x) is the same 
for each vector multiplet. To leading order in N, the matrix model for the A/" = 3 necklace 
theories involves extremizing a free energy functional of the type 



F[p, 5ya] = TiN'/^ / pix)dx 



UfIx] + 2x ^ QaSVaix) + p{x) ^ f{5ya{x)) 



a=X 



a=l 



(3.4) 
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where 5ya = Ua-i — Va, = J2a''^a, and / is a periodic function with period one given by 



fit) 



t when 



- < t < - 
2 - - 2 



(3.5) 



This free energy should be extremized over the set 



C = |(p,5i/a) : j dxp{x) = l;p(x) > and ^^p{x)6ya{x) = a.e.j 



(3.6) 



where we think of p{x) and p{x)6ya{x) as functions defined almost everywhere (a.e.). To 
enforce these constraints, we introduce the Lagrange multipliers p and z/(x): 



F[p,Sya] = F[p,Sya]-2nN^/' 



/i dx p{x) — 1^ + J dx p{x)u{x) 5ya 



X] 



The equations of motion that follow from this action are 

d 



1 



X] + i^ciX - v{x)) 5ya{x)] = fi- -nF\x\ 

a=l 

\f{^ya{.x))p{x) + qaX = U{x) . 



(3.7) 



(3.8a) 
(3.8b) 



The solution of these equations and the constraint Y^a=i ^Va = is 



p(x) = sl{x) - ss{x) , v{x) = - - [sl{x) + ss{x)] 



1 Isiix) + qax\ - \ssix) + qax\ 

2 sl(x) - ss(x) 



where we denoted by sl{x) and ss{x) the two solutions of the equation 



-nplxl + o 2^ \s{x) + QaX] = p , 



a=l 



(3.9) 



(3.10) 



with sl{x) > ss{x). Note that the set of s and x satisfying (3.10) defines a polygon. In fact, 



it defines the polygon V of (2.21) but rescaled by a factor of 2p. 



We show in Appendix [B] quite generally that in the continuum limit the extremized value 
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of the free energy F is proportional to the Lagrange multipler /x: 



4vriV3/2 



-/X, 



(3.11) 



Thus to determine F, it suffices to find /i. 



By definition, the density p(x) should integrate to one. The solution (3.9) demonstrates 



that the density p{x) is proportional to the length of a slice through V at constant x. Thus 
integrating p{x) over x should yield a quantity proportional to the area of the polygon. We 
obtain 



dx p{x) = 4/^2 Vol(P) 



(3.12) 



Assembling ( |3.1l| and yields F = 27rN^/^ /3^/Vo\{V). From Corollary |2| and in 



particular (2.25), we recover our AdS/CFT prediction (1.3). 



3.2 Operator Counting and the Matrix Model 

In this section, we relate the matrix model quantities p{x) and p{x)6ya{x) to numbers of 
operators in the chiral ring that don't vanish on the geometric branch of the moduli space 
of the A/" = 3 necklace quiver in the abelian case = 1. To that end, we first characterize 
the chiral ring, defined to be gauge-invariant combinations of the bifundamental fields and 
monopole operators modulo superpotential relations. At arbitrary A^, the superpotential for 
these M = 3 necklace quivers when ni? = is 



d 

W = J2v Tr(5„+iA„+i - AM 

^ 1 " 

a=l 



(3.13) 



where ka = qa+i — qa- When A^ = 1, this superpotential gives rise to d—2 linearly independent 
relations on the geometric branch of the moduli space 



(3.14) 



(As pointed out by |20|, these relations are one set of moment map constraints in the hy- 
perkahler quotient discussed previously.) There also exists a monopole operator T and 
anti-monopole operator T that create one and minus one units of fiux, respectively, through 
each gauge group. These monopole operators satisfy the quantum relation TT = 1 23,24 
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The Aa and Ba fields have i?-charge 1/2, guaranteeing that the i?-charge of the superpoten- 
tial is two. We will take the monopole operators to have zero /2-charge. Gauge invariance 
for this U{lY gauge theory means that the total U{lY charge of a gauge-invariant operator 
constructed from Aa, Ba, and monopole fields will vanish. Let Aa (Ba) have charge +1(— 1) 
under gauge group a and charge — 1(+1) under gauge group a — 1. The monopole operators 
have gauge charges ±ka under the ath gauge group. To write a gauge-invariant operator in 
a compact form, it is convenient to define the operators 

C^mg,+s ifmg„ + s>0 Tt™ ifx>0 

^rnqa+s = } ■ = { . (3.15) 

^^-mqa-s iimqa + s<0 [T"" If X < 

The gauge-invariant operators are 

0{m, S, I,]) = Jjm^^1l+s^n^g,+s _ _ _ ^^^^^ )i _ ^g^^g) 

Note we have used the superpotential relations to eliminate all but two of the (AiBi) factors 
that could potentially appear in a gauge-invariant operator. 

Let us consider the set of operators of the form 0{m, s, 0,0) with i?-charge less than r. 
The number of these operators is equal to the number of lattice points inside the polygon 

Vr = l{m,s) eR^ -.-^{s + qam\ <r\ . (3.17) 



a=l 



In the large r limit, the number of these lattice points is well approximated by Area(Pr)- The 
polygon in Section 2.2 and this polygon are related via V1/2 = V. Corollary [2] established 



that Area(P) and Vol(y) are proportional. As Area('Pr)=4T'^ Area(P), we have an additional 
relation between the number of a certain type of operator and Vol(y). 

The relation between Area('P) and an operator counting problem reveals an additional 
relationship between the eigenvalue density p{x) and the chiral ring. We claim that in the 
large r limit the number of operators 0{m, s, 0, 0) of i?-charge less than r and monopole 
charge between m = xr/fi and m + dm = (x + dx)r/fi is 



^2 

—p(x)dx . (3.18) 
1^^ 



From (3.9), the quantity p{x)dx in the matrix model corresponds to the area of a constant x 



strip of Vfj, of height sl{x) — ss{x) and width dx. Roughly speaking, the operator counting 



19 



gives p{x) a new interpretation as the number of operators 0{x, s, 0, 0) of i?-charge less than 
fi and monopole charge bounded between x and x + dx. As /i is of order one in the matrix 
model, we should be more careful and consider first Vr for some large r, giving rise to the 
factors of r//i in the claim. 

There is a subtle relation between 6ya and the chiral ring where Aa or Ba is set to zero. 
In the large r limit, we claim that 



^ (^SVaix) + ^ ) pix)dx 



(3.19) 



counts the number of operators of the form 0{m, s, 0, 0) with i?-charge less than r and with 
monopole charge between m = xr/p and m + dm = {x + dx)r/fi and with Ba = 0. Flipping 
the sign of 6ya{x) yields an equivalent expression for operators with Aa = 0. Note that 



Pix) Sya{x) + 



QaX < -Sl{x) , 

Sl{x) + QaX -Sl{x) < QaX < -Ss{x) 

Sl{x) - Ss{x) -Ss{x) < QaX . 



(3.20) 



For simplicity, let us first assume that /i is large and count the operators with i?-charge less 
than /i. We get the correct counting when g^x < —Sl{x) and when —Ss{x) < qaX because the 
region qaX + Ss{x) > corresponds to a portion of the polygon where the 0{x, s, 0, 0) contain 
no Ba while the region qaX + sl{x) < corresponds to operators that contain no Aa. In the 
central region, the operators that contain no Ba satisfy the constraint Sl{x) > s > —qaX. 
Because there is one operator per lattice point, the number of operators that contain no Ba 
is proportional to the difference sl{x) — {—qax). Similarly p{x){—5ya{x) + 1/2) will count 
the number of operators with no Aa- Again, since p is of order one in the matrix model, we 
should rescale our results for a polygon Vr in the large r limit, yielding the extra factors of 



r/ p in the claim (3.19) for 6ya- 

For a general supersymmetric gauge theory, we will not be able to make a clean separation 
between all operators in the chiral ring and operators of this special form (9(m, s,0,0). It 
is thus useful to reformulate these statements about p{x) and 5ya in terms of all chiral 
operators. Let us introduce the function ip{r, m) which counts the number of operators with 
i?-charge less than r and monopole charge less than m. We claim that 



dr'^dm 



m=rx/fi 



-p{x) 



(3.21) 
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Note that in the large r hmit, {d'^ip / drdm)dr dm can be interpreted as the number of opera- 
tors of -R-charge between r and r + dr and monopole charge between m and m + dm. Given 
an operator 0{m, s, 0, 0) of i?-charge ^ Yla \ ^ + Qa^l = i^q < r, we can form an operator of 
i?-charge equal to r by multiplying C(m, s, 0, 0) by a factor of {AiBiY {A2B2Y~^^''~^ . There 
are precisely r — tq + 1 ways of forming such a factor (assuming r — tq is an integer). This 
multiplicity r — tq + 1 associated to a lattice point (m, s) of the polygon can be interpreted 
as the total number of operators of fixed m and s with i?-charge equal to r. The difference 
in the number of operators with i?-charge r + 1 and r integrated over a strip at constant m 
now has the dual interpretation as rp{fim/r)/^ or d^ip / dr'^dm (in the large r limit). 

As for 5ya{x), let us introduce ipAai.f^''^) and ipBaij'i^) as the number of operators with 
i?-charge less than r and monopole charge less than m, with = and -Ba = respectively. 
We claim that 



drdm 



drdm 



m=rx/ II 



m=rx/ II 



%{x) (^-5ya{x) + ^ ) , (3.22) 

r ( 1 

-p(x) \ bya{x) + - 



''-pix) ( 5y„(x) + ^ ) . (3.23) 



Note that [d'^ipBal 9rdm)dr dm can be interpreted as the number of operators of -R-charge 
between r and r -|- dr, monopole charge between m and m + dm, and no Ba operators, in the 
large r limit. Given an operator 0{m, s, 0, 0) of -R-charge tq < r that does not involve the 
field Ba, we can form a unique operator of -R-charge equal to r by multiplying 0{m, s, 0, 0) 



by {AbBhY''^'' where h ^ a (provided r — tq is an integer). Thus, eq. (3.19) is counting the 
number of operators in the chiral ring with -R-charge between r and r -|- dr, magnetic charge 
between m = rx/p and m + dm = r{x + dx)/fi, and no Ba operators. There is a parallel 
argument for ipAa- 

3.3 Operator Counting and Volumes 

Given the close relation between p{x) and Vol(y), it is not surprising that there is also a close 
connection between p{x) and numbers of operators in the chiral ring. Motivated by Weyl's 



Law for eigenfunctions of a Laplacian on a curved manifold, the authors of 34 noticed that 
the number of holomorphic functions on a certain class of Calabi-Yau cones could be related 
to the volume of the Sasaki-Einstein manifold base. Their result holds for a Calabi-Yau cone 
which is a C*-fibration over a variety in weighted projective space with a Kahler-Einstein 
metric. This special case was later generalized to any complex Kahler cone by Martelli, 
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Sparks, and Yau 35 . Note that the set of holomorphic functions on the cone is precisely 
the chiral ring of the = 1 gauge theory and that the ring has a natural grading from the 
i?-charge. The relation between the number of holomorphic functions and the Vol(y) is 



Voi(r) 



1 



lim — lim ■?/'(r, m) 



r— >oo T m— ^-oo 



(3.24) 



where dim(y) = 2n — 1. In fact the result (3.24) helps to explain not only the relation 



(3.21) between p{x) and i/j{r,m) but also the relations (3.22) and (3.23) between 6ya{x) and 
restrictions of the chiral ring to rings where Aa = or Ba = 0. The subspace = of 
the Calabi-Yau cone is also a complex cone with Kahler structure. The level surfaces of this 



complex cone will be Sasaki manifolds which satisfy (3.24) with n = 3. The identifications 



(3.21), (3.22), and (3.23) along with (3.24) imply that 



p{x)dx 



p{x) { -6ya{x) + ^] dx 



p{x) ( 5ya{x) + - ] dx 



^Voi(y.j 

7^Vol(F.J 



(3.25) 
(3.26) 
(3.27) 



The relation (3.25) we derived already, but the second two relations are seemingly new. From 



the free energy functional (3.4), it is clear that 6ya{x) is an odd function of x, and thus for 



these necklace quivers the integral J p{x)6ya{x)dx vanishes trivially. The remaining integral 
yields the result 



Voi(y4j = VoKy^j 



TT 



Voi(r) 



(3.28) 



a result Yee found by other means [33] 



A field theoretic interpretation of these five dimensional cycles was provided by 41 , 42 



In an AdS4^ x Y solution of eleven-dimensional supergravity, an M5-brane wrapping such a 
cycle in Y looks like a point particle in AdS^ with a mass proportional to the volume of the 
cycle times the tension of the five-brane. The AdS/CFT dictionary provides a relationship 
between the mass of the particle and its conformal dimension. As the wrapped five-brane is 
supersymmetric, the conformal dimension can be related to the i?-charge of the corresponding 
operator that creates the state. If in the geometry the five-cycle corresponds to setting 
Aa = 0, the corresponding baryonic-like operator should involve an anti-symmetric product 
of A^ copies Aa- For our purposes, the essential point is a relation between VoI(Ya^) and the 
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i?-charge of Aa- 



6 Vol(F) 



(3.29) 



This point suggests a way of generalizing (3.22) and (3.23) to an arbitrary quiver gauge 



theory. We should replace the 1/2 with the i?-charge of the corresponding bifundamental 
field Xab with charge +1 under gauge group b and charge —1 under gauge group a: 



drdm 



m=rx/ fi 



p{x) [yb{x) - ya{x) + R{Xab)] 



(3.30) 



4 The Matrix Model for [p^ (/)-Branes 



As discussed in section 2.3, the free energy for a system of d (pa^Qa) five-branes can be 
computed on the gravity side by combining the M-theory prediction (1.3) with the eq. (2.22), 
where n = 2 and {pa, qa) = 0^- More explicitly, the polygon V takes the form 



V 



(x, s) e 



d 
a=l 



\PaS + gaa;| < 1 



(4.1) 



and the volume of V is related to F through F = 27rA^^/^/3 a/Vo1(P) . In the previous section 
we were able to reproduce this formula from the field theory side in the cases where Pa = 
or Pa = 1 for each a, as these were the cases where a simple Lagrangian description of the 
field theory led [T] to an expression for F in terms of a matrix integral. For > 1, as 
mentioned in the introduction, a Lagrangian description would involve the coupling to the 



T{U{N)) theory described in 30 . We will now use our large N intuition to figure out the 
matrix model at finite A^. 

As a first step towards obtaining such a matrix model, we note that at large A^ one can 
obtain the correct value for F by extremizing the free energy functional: 



F[p,6ya] = 7iN''/^ / p{x)dx 



2X —Sya{x) + p{x) Paf{Sya{x)/Pa) 

a=l a=l 



(4.2) 



The saddlepoint is a generalization of the earher eigenvalue distribution (3.9) for the (1, g^)- 
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-1+1/k 



1/k / i-i/k'^ 



a) b) c) 

Figure 3: a) The polygon for ABJM theory which can be built from a (1, 0) and (1, k) brane; 

b) the S-dual configuration involving a (0, 1) and {—k, 1) brane; c) an SL{2, Z) transform to 
a (1, 1) and (1 — k, 1) brane. 




branes: 



p{x) 



Sl[X) - Ss[X) , 
1 \PaSLix) + qaX\ - \PaSsix) + QaX 



u{x) = [sl{x) + ss{x)] 



(4.3) 



sl{x) - ss{x) 



where sl{x) > ss{x) are the two solutions of ^ Yl'!i=i \Pa.^ + 1a.x\ = Z^- 



The equation (4.2) seems ill behaved in the limit Pa — > 0, but in fact it is not. What 
happens is that 6ya{x) will saturate very quickly to ^Paf^ as we move away from a; = 0. 
Thus the function f{Sya/pa) will vanish, and 2xqaSya/pa can be replaced with qa\x\. Before, 
we identified the number of flavors with the sum np = Yip =o 9'"' ^^^^ 



a term of the form np\x\ will appear in the free energy (4.2). 



To give the reader a better sense of how SL{2, Z) acts on the polygon, we show V for 
ABJM theory and some of its SL(2, Z) transforms in figure [3] While all three polygons have 
the same volume (as they must given Corollary [3]) , the eigenvalue distributions can look 
quite different. In figure [3^, p{x) is a constant function for — < x < l/Zc. In figure |3]3, 
p{x) has two piecewise linear regions for —1 < a; < 1. Finally, in figure p{x) has one 
constant region and two linear regions in the interval — 1 + 1/A;<x<l — l/Zc. 



4.1 gj-branes at finite N 



One can go a little further and conjecture a finite N analog of the matrix integral (3.1 ) that 



reduces to (4.2) in the large limit. If we move the sinh and A^! factors from (3.2) to (3.3) 
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then the contribution from a (l,ga)-brane to the partition function is given by 



exp 



TTiga " ^1* ^ 



X 



Yl-^. 2 sinh 7r(Aa_i,i - A^-ij) ni<,- 2 sinh 7r(Aa,i - A^, 



(4.4) 



nij2cosh7r(Aa-i,i - A^,, 
We argue that the generahzation to {pa, qa) branes is given by 



L{pa,qa){^a-l, Aa) 



\Pa\ 



-N 



exp 



{y1 ^'-1.* - Yl ^li^ 



U^<, 2 Sinh ^(A,_i, - A,„i,,) n.<, 2 sinh ^(A,, - A,,,) 



(4.5) 



Pa 



<1 



Pa 



ni,,2cosh^(A._i,-A 



where this formula is correct only when {pa, Qa) are relatively prime. An {npa, nqa) five-brane, 
with (pa, qa) relatively prime, should be thought of as n (pa, qa) five-branes. Eq. ( |4.5 ) is based 



largely on the structure of (4.2). The term involving the hyperbolic cosine in (4.4) gives rise 



to the second term in (4.2) in the case Pa = 1. As the imaginary parts of the eigenvalues 



appear with a factor of pa in (4.2), at finite A^, we should divide the eigenvalue differences 
Xa-i,i — \a,j by Pa- The exponential term in (4.4) gives rise to the first term in (4.2) in 
the case Pa = 1. As a result at finite A^ we should replace the coefficient qa with qa/Pa in 



(4.4). This replacement suggests that the field theory on (p^, ga)-branes in some sense can 



be thought of as having Chern-Simons couplings -^{qa+i/pa+i — qa/Pa) J tr A dAa- (Such 

Recall that the factor exp[7ii{qa+i — qa) J2i il in 



an observation was also made in 
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the original matrix model comes from the classical contribution of the Chern-Simons term 
4^ ila+i — qa) f tr A dAa and its supersymmetric completion. 



The remaining factors of pa are required for (4.5) to be invariant under S-duality as we 



will now see. As an added bonus, by studying the action of S-duality, we will be able to 
deduce the partition function for the T{U{N)) theory of |30| . 

Our {p, g)-brane construction exists in type IIB string theory which is well known to be 
invariant under the action of SL{2, Z). One of the generators of SL{2, Z) is S-duality which 
we define to map a (p, g)-brane to a (— g,p)-brane. The work of Gaiotto and Witten pO] 
suggests that we also should be able to realize S-duality locally, on one (p, g)-brane at a time. 
The fact that S-duality squares to minus one suggests that we may be able to realize it as a 



Fourier transform acting on (4.5) 



For simplicity, we will restrict to the case where the ranks of the gauge groups are equal 
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to A^. Before introducing the Fourier transform, we make use of the identity 

Ui<j sinh(xj - Xj) Yl-^j sinh(|/j - yj) 
Ui,j cosh(xi - Vj) 

^ sech(a;i — yi) sech(x2 — Vi) • • • sech(a:„ — yi) ^ 

sech(xi - 2/2) sech(x2 - 2/2) • • • sech(a;„ - ^2) 

= det 

: : dots 

y sech(xi — ?/Ar) sech(x2 — i/Ar) ... sech{xN — yN) J 



(4.6) 



previously considered in 43 . Given this identity, we can write (4.5) as a sum over permu- 
tations 



/ \2 2 \ 



1 



2cosh|(Aj -cTpQ)) 



(4.7) 



We claim that a local S-duality is implemented by the following Fourier transform: 



To demonstrate this claim, we isolate the integrals over \j and crp(j): 



\P\J ' cosh^(A,-a,(,)) • 

With the change of variables x± = \j ± <7p{j), this integral is straightforward to perform: 



(4.9) 



2\p\ 
1 



dx-L dx- e 



Tlx. 



P 



, , , c/x_(5(^x_+/i, -z/„(,) )e™-('^^+^''0))sech — 
bl i VP 7 P 



(4.10) 



exp 



•P/ 2 2 ^ 



1 



coshf(/ij - z/p(j)) 



Taking the product over the eigenvalues and averaging over permutations yields (4.8). 



This local S-duality composes in a nice way. Consider applying similar Fourier transforms 
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to neighboring (p, gf)-branes: 

= y e2-^-^L(,,,)(A,/i)L(,,,,)(/.,i>)e-2-'^-d^Ad^z>d^^ (4.11) 

Thus if we apply a local S-duality to each (pa, ga)-brane in the necklace, the factors of e^'^'^'^'^ 
cancel out and the resulting partition function is invariant under a global action of S-duality. 

We would like to give a better interpretation of this group action. Consider acting on a 
single (p, g)-brane with this local S-duality: 

j L(_,,^)(A,/.)L(^,,,,)(/.,i/)d^/. = j e2'^^^-^L(p,,)(A,/i)e-2"'^-'^L(^,,,,)(/.,i/)(i^Ad^/id^/x. 

(4.12) 

One way of interpreting the e~'^'^^^'^ is to posit that some new object has been inserted 
between the {p, q')-brane and the {p', g')-brane that implements a local S-duality. This object 
contributes to the partition function an amount 

L5(//,/i) = e-2-''-^. (4.13) 

Similarly, right before the {p, q) brane we introduced another object that undoes the local 
S-duality: 

L5-i(A,A) = e2'^^^-\ (4.14) 

Let us see how the T{U{N)) theory arises. So far, we have been thinking of the ip,q)- 
branes as the building blocks out of which we construct the partition function. Alternately, 
we can decompose the partition function into the contributions from the D3-brane segments 
and associated U (N) gauge groups. Prom the D3-brane point of view, the object I/5(//, fl) 
implementing S-duality splits a D3-brane segment into two regions, each characterized by a 
U{N) gauge theory, one with eigenvalues /i and one with eigenvalues jl. The -£'(p,g) and -£'(^',5') 
factors do not have enough factors of hyperbohc sine to describe two U (N) gauge theories. 
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These factors are simple to recover if we say the new object has a partition function 



This object has a natural interpretation as the partition function of the T{U{N)) theory]^ 



This partition function has been found independently by 31 . 

At this point, it is clear that we should be able to realize any element of SL{2, Z) acting 
locally on our necklace theories. The group SL{2, Z) has two generators: S, which we 
discussed above, and T. We define T to send q ^ q + p and leave p invariant. If we 
think of {p, q) as a two-component column vector on which SL{2, Z) acts in the fundamental 
representation, then S and T are the two-by-two matrices 





They satisfy S"^ = —1 and {ST)^ = 1. To find the action of T on the matrix model, one can 



see from (4.7) that 

L(,,,+,)(A,a) = e'^*^-^-^L(,,,)(A,o-)e-™-'^+^^ , (4.17) 

where 6' is a phase to be determined. Similarly to Ls, one can therefore define 

Lr(a) = e-"*"-"+*^ , L^-i (A) = e"*^'^"^^ , (4. 18) 

so one can describe the contribution of a {p,q + p) five-brane to the matrix model integrand 
as L'r-i(A)L(p^g)(A, cr)LT(cr), where the factor of Lt{<j) corresponds to the local action of a 
T-transformation and Lt-i{X) corresponds to the local action of T^^. 

We have defined the effect of the 5* and T generators on the matrix model so that S 
requires us to have two distinct sets of eigenvalues /i and ft in the two regions separated 
by an "S-boundary" , while T doesn't, acting just by multiplication of the eigenvalues a in 
the region containing the "T-boundary" by Lt{(t). We could have said, however, that a 
T-boundary also requires two distinct sets of eigenvalues cr and a around it, in which case 
we should have described its contribution as LT{(T)5{a — a). By the same logic, it follows 

^Actually, the T{U{N)) partition function needs to be antisymmctrizcd with respect to permuting /i. But 
since L(p g)(/x, A) is aheady antisymmetric in /x, we can get away with not antisymmctrizing Lg. 
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that the identity element should be also viewed as a Dirac delta function Li{X, a) = S{X — a). 
The operator equal to minus the identity is L_i(A, a) = 6{X + a). We can verify by explicit 
computation that 



^ss- 



.(/i, A) 



2niaiX-^)^N^ = - A) = A) 



(4.19) 



and similarly that L^t-^ = Li, L52 = L_i and L 



iST)3 



Li provided 36* = -kN/A. 



We should mention that the expressions (4.13), (4.14), and (4.18) can also be justified 
by analyzing the case = 1, where the action of S'L(2,Z) on Chern-Simons theories was 



described in 30 ,44 . Indeed, as explained in 44 , the T generator just shifts the CS level 



by one unit, so the action changes by — A dA plus its supersymmetric completion. 

The classical contribution to the partition function from the scalar A in the M = 2 vector 
multiplet then gives Lt{X). The action of 5* on a CS theory with some gauge field A consists 
of introducing another dynamical gauge field A that couples to the topological current *dA. 
This coupling takes the form of an off-diagonal Chern-Simons term J A A dA plus its 
supersymmetric completion. If /i and fl are the scalars in the corresponding vector multiplets, 
the classical contribution to the partition function from this off-diagonal Chern-Simons term 



is precisely given by (4.13). 



5 Discussion 



Our main results are additional evidence presented in Section 2.3 for the F-theorem conjec- 



ture, the relations (1.4) and (1.5) between numbers of chiral operators and eigenvalue distri- 



butions, and a conjectured form (4.5) of the matrix model corresponding to a {pa, ga)-brane 
construction in type IIB string theory. Each of these results requires some brief discussion. 



We would like to investigate further our proposed matrix model (4.5). In particular, it 



would be interesting to see how (4.5) transforms under Seiberg duality. One statement of 



the s-rule 45 is that a theory breaks supersymmetry for which there exists a Seiberg duality 



that produces a gauge group with a negative rank. Given the matrix model's status as a type 
of supersymmetric index, one expects that the partition function should vanish for theories 



that violate the s-rule 46 



Regarding the F-theorem, we have not constructed any explicit RG flows, either on the 
gauge theory side or, via the AdS/CFT correspondence, on the gravity side. Instead, we 
have posited the existence of reasonable seeming flows, and we have examined F at the IR 
and UV fixed points. For example, by adding a mass to fundamental flavors, one should be 
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able to flow from a theory with a (0, g)-brane in the UV to one without it in the JR. The 
corresponding volume of the tri-Sasaki Einstein manifold will increase, leading to a decrease 
in F. Similarly, we can consider an RG flow where a (p, g)-brane forms a bound state with 



a (p', g')-brane. Under such a flow, F will also decrease. Given the result of refs. 15 ,16 



described in the introduction, it seems likely that any gravity dual of an RG flow will obey 
the F-theorem. One way of interpreting our results, given that our flows also obey the F- 
theorem, is that it may be possible to realize these flows as solutions of eleven-dimensional 
supergravity. 



Given that the operator counting relations (1.4) and (1.5) can be defined for essentially 
any KWY matrix model and corresponding superconformal field theory, one wonders if they 
hold more generally. In a sequel to this paper [47j, we investigate the large limit of 
the KWY matrix models for a number of other superconformal field theories, and we find 
that indeed these relations are always satisfied. We look at necklace quivers with additional 
adjoint and fundamental fields. We look at a couple of non-necklace quivers, for example a 
Z2 X Z2 orbifold theory. In order to look at theories with M = 2 supersymmetry, we have to 
generalize the KWY matrix model to allow for arbitrary /2-charges [2]. It turns out that in 
the theories we study, relations (1.4) and (1.5) are valid not just for the correct i?-charges 
but for any i?-charges compatible with the marginality of the superpotential. 

One constraint in these investigations is that for chiral theories, the KWY matrix model 
does not seem to have a large limit that is compatible with a dual eleven-dimensional 
supergravity description [5]. It will be interesting to see if the relations (1.4) and (1.5) can 
give any insight into how the matrix model might be modified to allow for such a limit. 
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A Proof of the tree formula 



We prove the Tree Formula by assuming that the columns are ordered in the way described 
in Corollary [3j Thus, the proof reduces to showing that the following equation holds 



E 



7a(a+l) 



We have the "Ptolemy relations" 

gabled = laclbd + ladlbc if ab / / cd. 

The notation / / means a and h separate c and d. We can use the relations to show 

l(a-l)(a+l)(^a + 27(a_l)a7a{a+l) = l(a-l)a(^a+l + 7a(a+l)0"a-l- 



(A.i: 



(A.2) 



(A.3) 



Our starting point is the Kirchhoff matrix-tree theorem which gives a relation between 
the absolute value of a certain determinant and the sum over trees. In particular, consider 
the matrix Q* where 



0"2 


-723 


-724 • • • 


-l2d 


-732 


0-3 


-734 • • • 


-l3d 


-742 


-743 


0-4 


-lAd 






: dots 




-ld2 


-7<i3 


~ldA ■ ■ ■ 





(A.4) 



The matrix-tree theorem states that 



detQ*= n 

(V,E)£T {a,b)eE 



(A.5) 



We observe that if we take 7(a+i)a times the (a — 2)nd row minus 7(a_i)(a+i) times the 
(a — l)st row plus 7a(a-i) times the ath row, then most of the entries will cancel out. So we 



31 



have AQ* = B where 



A 



( 7l3 + 723 —712 + 723 
-734 724 

-745 



723 
-723 
735 



723 






B 



\ l{d-i)d l{d-i)d l{d-i)d 

( 7l20"3 + 7230"! -7l20"3 

-734Cr2 7230"4 + 734O" 2 

-745(^3 



(A.6) 



dots : 

• • • 7(d-i)i + l{d-\)d j 

... \ 

-7230-4 ... 

7340-5 + 7450-3 • • • 

: dots : 

... i{d-\)d(y\ + id\(yd-\ ) 



A.7) 



In constructing A we used the fact that the missing row of Q is minus the sum of all of the 



other rows. We used (A. 2) and (A. 3) to simplify B. Tri-diagonal matrices satisfy a three 



term recurrence relation. If we let B^ be the matrix consisting of the first a rows of -B, then 

det Ba = (7a(a+l)0-a+2 + 7(a+l)(a+2)0-a) dct - 'y(a+l){a+2}l{a-l)a<^aO-a+l det Ba-2- (A. 8) 



We can show by induction that 

det{Ba) = 
and in particular 

det (5) = det(Erf_i 



a+2 \ / a+1 



(A.9) 



b=2 



v6=l 



7l27dl 



n ^a7a(a+l) Yl 



,a=l 



la{a+l) 



(A.IO) 



Therefore we want to show that det (A) = 2cri 11^=2 7a(a+i)- 

We let A' be A with the 723 in the top row and the 'y(d-i)d in the bottom row removed. 
Let A'^fj be the matrix containing the ath through 6th rows and columns of A'. Write the 
determinant of A as 

det (A) = {A[ + 723M) A A'2 A . . . A A'a_^ A + 7(d-i)dM) 



where u = (1,1 



, . . . , 



1) and A'^ are the rows of A'. Expanding out the anti-symmetric product, 
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we find that 



det{A) = det(A') + 7231^ A ^'2 A ■ ■ ■ A K-i + 7(d-i)d^i A ... A A'a_^ A u 



= det(^') + E 



a=l 
d 



ia-\ 



723 



n7(6+2)(6+3) det^J„+i)(^_ 



1) 



+7<i(d-i) ( n 7(6-i)(6-2) ) det 

\h=a+2 



(A.11) 



Using the three term recurrence relation for the determinant of a tridiagonal matrix, we 
can show by induction that 



det = 7i(a+2) n ^^^^^ 



1) ' 



6=2 



d-2 



det^a(d_l) = 7laJJ 7(6+1) (6+2) , 



(A.12) 



b=a 



det A' = 0, 



(A. 13) 



where last equation follows from the second one by setting a — 1 and using 711 = 0. Therefore 

det(A)=2(E7iJ (n^»(»+i) 

\ a / \a=2 

d-l 

= 2ai Yl 7a(a+l) 
a=2 

as expected. 

B A virial theorem for matrix models 

To leading order in N, the matrix model for superconformal field theories with M-theory 
duals involves extremizing a free energy functional of the type 

F[p,y„] = j dxp{xff{ya{x)) - j dx p{x)V{x,ya{x)) , (B.l) 
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for some functions / and V. In all examples, V is homogeneous of degree 1 in x, namely it 
satisfies 



V = x 



dV 
dx 



(B.2) 



The free energy functional (B.l ) should be extremized under the constraint that p is a density 



normalized so that J dx p{x) = 1. Such a constraint can be implemented with a Lagrange 
multiplier /x by defining a new functional 



F[p,y,] = F[p,ya] - 27iN^/'p (^j dx p{x) - 1^ 



(B.3) 



and varying ( B.3[ ) instead of ( B.l[ ) with respect to p and ya- If we denote by F the on-shell 
value of F[p, ya] (or of F[p,ya]), we will now show that 



F = ^—p, 



(B.4) 



regardless of which particular matrix model we're solving provided that eq. (B.2) is obeyed. 



B.l A slick proof 



The equations of motion following from (B.3) are 



2p{x)f{ya{x)) - V{x,ya{x)) = 27cN^/'p., 
p{xfdaf{ya{x)) - p{x)daV{x,ya{x)) = 0, 



(B.5) 



Differentiating the first equation with respect to x, multiplying it by and using the 

second equation we get 



p{x) 



2p'{x)f{ya{x)) + p{xf J2 daf{ya{x))y'j 



X] 



d 

p{x) — V{x,ya{x)) . (B.6) 



To prove (B.4), consider the function 



G{x) = xp{xff{ya{x)) 



(B.7) 
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Let's calculate the derivative of this function with respect to x\ 



G\x) = p{xff{ya{x)) + 2xp{x)p\x)f{ya{x)) + xp{xf ^ daf{ya{,x))y^{x) 



d 

P{xff{ya{x)) + Xp{x) — V{x, ya{x)) 
P{xff{ya{x)) + p{x)V{x, Vaix)) , 



(B.8) 



where in the second line we used (B.6) and in the third line we used (B.2). The function 



G{x) vanishes at ±00, so J dxG'{x) = 0, which from (B.8) implies that on-shell we have 

dx p{xYf{ya{x)) = ~ I dx p{x)V{x,ya{x)) . (B.9) 



From (B.l) we have then that 



F = 2 dxp{xff{ya{x)). 



Multiplying the first equation in (B.5) by p[x) and integrating in x we get 



27rAr3/V = 3 / dxp{xff{ya{x)). 



Taking the ratio of the last two equations one obtains (B.4). 



(B.IO) 



(B.ii: 



B.2 A more enlightening proof 

The proof given above is really that of a virial theorem. To put the virial theorem in a 
more familiar form, let's define the cumulative distribution t[x) = dx'p{x') G [0, 1] and 



express (B.l) as the action 

S[x{t),ya{t)] 



dt L , 



L 



fiVa) 



X 



V{x,ya) 



(B.12) 



where x and ya should be thought of as functions of t. To go to a Hamiltonian formulation, 
we introduce the momentum conjugate to x: 



Px = 



dL 
dx 



fiVa) 



(B.13) 



x^ 
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The Hamiltonian is 

H = p,x-L = KE + PE, KE = = -2f{y,y/',/^, , PE = V{x, y^) . 

X 

(B.14) 

In deriving a virial theorem for this Hamiltonian, one can just ignore the dependence on ya 
because the ya are non-dynamical. Since the kinetic energy KE is homogeneous of degree 
1/2 in px and the potential energy is homogeneous of degree 1 in x, we have 



-(KE) = (PE) 



(B.15) 



where (■ ■ ■ ) means J^dt {■■■). From ( |b12| one sees that 



F = S-on-sheii = -^{KE) - (PE) = -(KE) . 



From the first equation in (B.5) we see that H = —2nN^/'^fj., so 



— = -{H) = -(KE) - (PE) = -|(KE) 



Taking the ratio of (B.16) and (B.17) one again obtains (B.4). 



(B.16) 



(B.17) 
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